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ABSTRACT 


The purpose of this research is to design a digital phased array for receiving and 
tracking a 2.4 GHz frequency modulation (FM) video signal from an unmanned air 
vehicle. The tracking is done using a monopulse technique. Various numbers of elements 
were simulated to access the pattern coverage. The beamforming and demodulation are 
performed digitally by a computer. Several antenna architectures were studied and the 
requirements for hardware components such as the analog-to-digital converters and 
amplifiers are specified. Several methods for generating the carrier signal for the local 
oscillator are also discussed. They include an independent local oscillator, a voltage 
controlled oscillator with feedback (1.e., a phased locked loop) and mixing with a phase 


shifted copy of the received signal (Bilotti’s method). 
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EXECUTIVE SUMMARY 


As a part of an ongoing project, this research aimed at designing a digital phased 
array for receiving and tracking a 2.4 GHz FM video signal from an unmanned air 
vehicle (UAV). Tracking the signal from the UAV will establish a data link for data 
transfer such as video, image or audio between the UAV and the ground station. The 


tracking is done using a monopulse technique. 


An atray tradeoff study was done to determine the benefits of 1 row and 2 row 
arrays, with 8 elements per row, in terms of the half-power beamwidths (HPBWs) of the 
azimuth and elevation patterns of the array. The one row array was chosen primarily due 
to its wide elevation pattern coverage. This eliminates the need for elevation tracking as 
the UAV is within the HPBW for all UAV ranges of interest. Adding a row of 8 more 
elements to increase the signal-to-noise ratio (SNR) will increase the cost of equipment. 
The additional SNR is not required, as the UAV transmitter power is sufficiently high. 
The sum and difference beams for the Matlab simulation of the one row array are verified 
with the CST Microwave Studio simulations. The printed circuit dipole element was 
found to have an E-plane pattern HPBW between that of an ideal dipole and a half-wave 


dipole. 


Several approaches for the demodulation of FM signals are explored. The 
considerations for the choice of the approaches for the FM demodulation are the 
minimum circuitry requirement, the ability to conduct beamforming, the recovery of the 
transmitted signal and the susceptibility to phase errors and local oscillator (LO) drift. 
Three FM demodulation techniques are considered: Bilotti’s quadrature demodulator, a 
local oscillator with no feedback and a digital phase locked loop (PLL) quadrature 
detector. All three use digital beamforming and signal processing to demodulate the FM. 
Simulink models of the three FM demodulation techniques show satisfactory 


performance for the recovery of a two-tone signal. 


The simulations were done at a carrier frequency of 50 kHz instead of 2.4 GHz, to 


greatly save on the processing time. A LO with no feedback was chosen, as it uses the 


Xlil 


least circuitry with satisfactory performance for the recovery of the two-tone signal. 
Beamforming can also be performed easily as the raw in phase (I) and quadrature (Q) 
data are preserved. This technique also eases the computer processing. Since no error 
correction needs to be computed for the LO as it does for a PLL. Simulations were also 
done with relative phase errors introduced in the channels and a frequency offset in the 
LO without feedback. The Simulink simulations indicate that the recovery of the two- 


tone signal can still be achieved with the expected errors. 


Finally, the possible hardware components are briefly introduced and the 
minimum power for the receiver module was used to determine the maximum operating 


range of the array using the one-way Friis equation. 
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I. INTRODUCTION 


A. BACKGROUND 


Phased array development began around 1958, when early radio transmitters and 
radars used multiple radiating elements to achieve the desired antenna radiation patterns 
[1]. As technologies advance, existing systems are improved and upgraded. Sometimes, 
new systems are created when the existing systems become obsolete. Unmanned systems 
and array antennas, two technologies that improve continuously, have brought many new 
capabilities to military operations. The superiorities of array antennas are the digital and 
computerized processing capabilities that can be built into them. Digital phase shifting to 
form a beam, the electronic and rapid beam scanning ability, and its broadband potential 


are the advantages of array antennas. 


Since the Gulf War in the 1990s, the military has a new awareness of the 
importance of technological advances such as stealth ability, unmanned systems and 
global positioning systems (GPS) in improving the success of a war. One of these 
technologies is the digital phased array. It has many advantages, such as scanning rapidly 
without physically moving. Multiple, simultaneous beams can be achieved by the means 
of digital-beam forming (DBF). Error correction and self-calibration can be achieved 
easily in the signal processing. DBF also allows pattern control to direct the nulls in the 
pattern in the direction of noise jammers, and also suppresses clutter, which is a cause of 
degradation in range. DBF also gives the opportunity to compensate for the failure of the 
elements in the array. Different data rates, depending on the operating environment, can 
be achieved with DBF. The array antennas can be used in many applications, such as 
communications, data-links, radar (search and track), and electronic warfare (EW). This 


shows perhaps the most important advantage of these antennas: the multi-functionality. 


Unmanned aerial vehicles (UAVs), one important type of unmanned system, play 
a very critical role in military operations. UAVs have a variety of uses, such as 
surveillance and reconnaissance, target recognition, battle damage assessment, EW, 


search and rescue, and traffic monitoring. All of these can be conducted without 
1 


endangering the life of a pilot, as they do for manned systems. Important roles for the 
UAVs are troop support and carrying out reconnaissance and surveillance missions, 
which require maintaining a data-link with troops in order to send any data collected, 


such as video images or audio, back to the command center. 


During operations, it is imperative to continuously maintain a data and control 
link with the operator. This requires the ground station antenna to track the UAV so the 
antenna beam is pointed properly. The purpose of this research is to design and build an 
array to angle-track a UAV and, eventually, to accomplish the data transfer from the 
UAV to the ground station. The antenna type on the UAV is not a subject of this research 
and it will be assumed to be a low-gain omni-directional antenna. What is of concern is 


the array system that will be designed to achieve the purpose that was just mentioned. 


Why is an array antenna selected for this application? The narrow beam of the 
antenna must be kept on the UAV as it moves and maneuvers. The beams of an array can 
track without physically moving the antenna. Digital array systems can be reconfigured 
and different beamwidths can be used, depending on the UAV path and the data rate 
desired for the type of operation. Array antenna systems can track and follow more than 
one UAV at a time (multiple tracking) [2]. Figure 1 shows a High Gain Advanced GPS 
Receiver (HAGR) phased array antenna. 





Figure 1. . HAGR 16 element antenna array. From [3]. 


B. PHASED ARRAY TECHNOLOGY FOR MILITARY APPLICATIONS 


Since 1901, array antenna technology has made tremendous progress with types 
that vary from linear to planar to conformal. Digital array applications enable the main 
beam of arrays to be steered electronically, rather than mechanically. A phased array has 
its radiating elements excited by the relatively phased currents and its main beam can 
scan the desired field continuously. Since the 1980s, phased arrays have gained wider use 


in radar and communication applications. 


Phased array radars are used in a number of navy ships, such as the US Navy 
Aegis ships, the German Navy’s F124 frigates and LCF frigates of the Royal Netherlands 
Navy [4]. The advantages of the phased array radars used in communication systems are 


summarized below: 


1. Electronically scanning the beam and its multiple element system 
increases the reliability of the radar system, and maintenance is simplified by 
replacing the faulty elements with new ones without a need to return to the 
shipyard. An antenna can have a significant number of its elements destroyed 
and still be able to function, albeit at a reduced capability. Electronic arrays 
have no moving parts, do not suffer from vibration and are inherently better 
stabilized against ship or aircraft motion, as the electronic array simply steers 


the main beam counter to the detected motion. 


2. It is suitable for low observable platforms as antenna elements of the array, 
like a microstrip patch antenna, can conform to the structure of stealthy 
platforms. Its agile scanning rate, with a narrow beam and lower sidelobes, 


gives a low probability of interception and is less vulnerable to jamming. 


3. The integrated circuit (IC) technology in microwave monolithic ICs 

(MMICs) is being driven by the need to integrate disparate technology 

functions into one unit. Through the use of IC technology, phased arrays can 

reach improved cost, size, and weight metrics that have previously prevented 

the widespread adoption in mass market scenarios. The high power handling 

capability and low loss of ferrite phase shifters allow power amplifiers, low 
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noise amplifiers and phase shifters to be chosen to achieve high processing 
performance, with relatively little impact to cost. Hence the solid state and IC 
technologies will benefit the phased array in terms of costs, weight and 


reliability [5, 6]. 


4. Its ability to create multifunction systems (radar, sensing, communications, 
etc.) with a single device enabled by the system’s software-driven nature 
increases the tactical flexibility of its operational employment. This capability 


is explained in the following scenarios: 


e The radiated power from each element can be adjusted when each 
element has its own transmitter. At the expense of cost, the controlling 
of the radiated power in sectors addresses the concern for emission 


control and co-site interference. 


e The elements can be programmed to communicate with multiple 
UAVs by scanning one single beam rapidly or producing multiple, 


independent beams. This would increase the dwell time on each UAV. 


e The multifunction system reduces the number of antennas on the 
military platforms, hence reducing weight and the radar cross section 


of the platform. 


e Its flexibility in electronic control allows for customized software 


solutions to address different user scenarios with one platform [5]. 


e A graceful degradation is allowed when some of the elements are 


damaged enabling the array system to continue functioning. 


The frequency bandwidth is a major design challenge for the array system as 
bandwidth is a scarce commodity in modern communications systems, with the Federal 
Communications Commission (FCC) deciding which parts of the spectrum can be used 
for emerging communication standards. Despite the multi-functionality of the current 
arrays, different functions need different bands, due to the unique characteristics and 


behaviors of the electromagnetic waves in different frequency bands. Multifunction 
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atrays can be used for surveillance, tracking, communications, direction finding and 
weapon guidance purposes. The radar should be designed in lower bands for longer range. 
The AN/SPY-1D operates in S band. Current electronic support equipment operates 
within the range from 2 to 20 GHz. However, due to the dense electronic environment, it 
is very difficult for the electronic warfare operators to detect radars operating in the UHF 


band in the littoral waters. 


The complexity of the phased arrays will increase the cost. Commercial off-the- 
shelf (COTS) hardware is to be exploited to reduce the cost of the phased arrays. The 
ultimate aim of the antenna engineer is to decrease the cost of the phased array antennas 


without sacrificing performance [2]. 


c. PREVIOUS WORK 


This research is part of the ongoing project that started with the thesis by Bartee 
[7]. Esswein [8] and Gezer [2] continued the work and demonstrated the feasibility of 


building a digital phased array with available COTS hardware. 
Bartee designed and constructed the distributed phased array at X-band. 


Esswein constructed a 24-element, three-dimensional, dual-ground plane, active 
phased array. A Genetic Algorithm (GA) was used to determine the element excitations 
and the measured radiation beam was in agreement with the theoretical calculations. The 
modulator feeding each printed antenna element used for this array was the Analog 
Devices AD8346EVAL quadrature modulator board, which is now being incorporated 


into a transmit/receive (T/R) module for radar applications. 


Gezer constructed a 1 by 8 phased array to track a signal transmitted continuously 
by a UAV operating at 2.4 GHz. This array used a monopulse difference beam to track 
the UAV signal. The array was designed based on typical UAV operational requirements, 
which led to the selection of the array parameters, such as the number of elements and 


spacing. 


D. SCOPE OF RESEARCH 


State of the art technology in electronic systems is a means to increase the 
information superiority and provide the decision makers with a well described tactical 
picture of the war scene. UAVs can carry out many reconnaissance and surveillance 
missions. They also provide the necessary information to the troops by being integrated 
into the network-centric warfare elements as seen in Figure 2. 
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Figure 2. The cooperation among the elements of a UAV system. From [10]. 


As an element of the network, the ground control stations of UAV systems can 
exploit the phased array technology. Array antennas are a very sophisticated means to 


accomplish accurate UAV tracking and maintain a good data-link. 


The research goal was to design a digital phased array for tracking and receiving 
FM video signals transmitted at 2.4 GHz from a UAV. The design incorporates COTS 
hardware to reduce cost. Simulation is another way to reduce the cost of the design and 
comprehend the consequences of design changes before manufacturing, and to see if 
there needs to be any adjustment. Thus, the design of the tracking array antenna system is 


verified through simulations before being constructed [2]. 
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E. ORGANIZATION OF THESIS 


Chapter II discusses the concepts of arrays. The concepts of a digital antenna and 
digital beamforming are explained, and the mathematical theory for arrays with ground 
planes is elaborated. Amplitude comparison monopulse systems, and how tracking of the 


UAV is achieved, along with external causes of errors, are mentioned. 


Chapter III addresses the results and analysis of the array tradeoff study. It 
discusses the results obtained for the 1 row and 2 row arrays with 8 elements. Also, the 


effect of polarization is studied. 


Chapter IV addresses the FM demodulation techniques being considered. Several 
methods for generating the carrier signal are also discussed. They include an independent 
local oscillator, a voltage controlled oscillator with feedback (i.e., a phased locked loop) 


and mixing with a phase shifted copy of the received signal (Bilotti’s method). 


Chapter V discusses the system architecture of the array and its considerations. 


Hardware implementation using COTS components is also described. 


Chapter VI summarizes what was done in this research and gives 


recommendations for future work. 
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Il. TRACKING ARRAY DESIGN 


A. WHY AN ARRAY ANTENNA? 


An antenna array is a configuration of individual radiating elements that are 
arranged in space to produce a directional radiation pattern. Single-element antennas have 
broad and low directivity radiation patterns, which are not suitable for long distance 
communication or accurate tracking. High directivity can be achieved by increasing the 
electrical dimensions of the single-element antennas in terms of wavelength, hence its 
physical size. This makes the antenna too cumbersome for mobile operations. The 
alternative is to combine multiple electrically small antennas to form an array. Antenna 
airays come in various geometrical configurations, including linear, planar and 
volumetric. The planar antenna array is often used in communication and _ radar 


applications, and therefore this general configuration will be analyzed [11]. 


Conventional phased arrays employ microwave circuits to distribute signals to the 
elements or, on receive, collect signals from the elements. These beamforming networks 
(or beamformers), as shown in Figure 3, can be bulky, narrowband and lossy. Some of 


these problems can be avoided by using digital beamforming. 
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Shifter 





Equal power 
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Figure 3. | Conventional microwave beamforming network. From [12]. 
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B. DIGITAL ARRAYS! 


The basic architecture of a digital array is shown in Figure 4. For a two-way 
communication application, each element must be capable of transmission and reception. 
In the transmit channel, the waveform is generated and amplified before passing to the 
array element. On receive, the signal is amplified and detected. The detected signals are 
collected by a processor which forms the beam and performs other communication 


functions. 


In order for all of the elements to send and collect data coherently, a common 
frequency (phase) reference must be distributed to all of the elements. This signal is 
referred to as the local oscillator (LO) signal. The elements must also be time 
synchronized so that transmission can occur simultaneously and the transmitted signals 


arrive at the other end of the link at the same time. 
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Figure 4. | Generic digital antenna architecture From [12]. 


Generally the array has a modular architecture where the transmit (Tx) and 


receive (Rx) functions are channeled. The transmit portion generates the waveform, 


! Much of the material in this section is adapted from unpublished notes by Prof. D. Jenn. 
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amplifies and filters it before passing it to the antenna. In modern arrays, the waveform is 
generated using a direct digital synthesizer (DDS). In its simplest form, a DDS is a clock 
and memory with a stored waveform. The waveform is read out at a rate determined by 
the clock. The DDS output is usually a baseband signal that has to be upconverted, 
although recent DDS clock speeds are achieving 2 GHz. Therefore, at VHF radar 
frequencies and below, it may be possible to amplify the DDS output directly without 
upconversion. For the tracking application, this research is concerned only with the 


receive function. Therefore, the remaining discussion only addresses the receive channel. 


On the receive side the returned signal is amplified with a low noise amplifier 
(LNA), downconverted, filtered and then sampled. State of the art analog-to-digital 
converters (ADCs) are in the 2 GHz range, so at VHF and below it may be possible to 


sample the waveform directly without downconverting. 


The in-phase (/) and quadrature (Q) form of the data are handled by the hardware. 


If the signal at carrier frequency @, =27 f, iS: 
s(t) = A(t) cos| 27 f,t + O(t)] (1) 


the in-phase and quadrature representation is: 


s(t) = I(t) cos(2z ft) - (0) sin(2z ft) (2) 
where 
I(t) = A(t) cos[@(z)] 3) 
O(t) = A@)sin[ O(n] 
and 
A(t) = 1(1)? + OW) 
(1) = tan”! gO. ~ 


The signal is conveniently shown in the //Q plane (Figure 5), where the distance from the 


origin is A(t) and the rotation from the J axis is O(f). 
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Figure 5 also shows a quadrature demodulator, the Analog Devices AD8347 [14], 
which operates from 800 MHz to 2.6 GHz with a 60 MHz instantaneous bandwidth. The 


demodulator has differential outputs (J seer gs Q*, and Q ) for obtaining the baseband 


in-phase and quadrature components of the received signal. The components are given 
by: 


T=l'-I- 


2 
Q=2°-Q 


The hardware aspects are discussed further in Chapter V, Section A. 


Q 


“be 
LY 


Figure 5. //Q signal diagram (left) and the AD8347 quadrature demodulator (right). 
From [12]. 
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The demodulation process is analyzed using the complex envelop of the signal 


u(t) [15]: 


u(t) =1(t)+ jO(t) 
s(t) = Re{u(nye/“*"| (6) 


1 ‘ . 
= sl uinel@et +u*(the 1c! 1: 
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The mixing process (shown in Figure 6) is equivalent to multiplying the received signal 


with 
i ae : 
COS O,.9 = 5 lei@ro! of ce J101 ) (7) 
in the J channel and with 
F 1 . . 
—sin@,9 = ~ 57 lei@ro! _¢ I10t) (8) 
J 


in the Q channel. 


LPF ~-—~» ADC ;-——/ 
r(t) 
90° 


a LPF ADC }—> @ 


Figure 6. | Quadrature demodulation process. From [12]. 









































For the in-phase channel: 


1 : : : : 
r(t) COS @LQ = qlueleet +u* (the Jee! ] (e/@Lot + e J10t) 


3 u(t) el (@t@poM 4 u(t) el(@c-@10)t 


: (9) 


ae *(t) e Jeo) 1 U oa e I(t @z0)t 


= Rey WO elie +@10)t 4+Re {Peltoe-eor| 


where Re{.} is the real operator and the * is complex conjugate. 


If direct down conversion is used (no intermediate frequency, @,, = @, ), then the 
first term is at frequency 2@,., which is filtered by the low pass filter (LPF). The second 


term is the desired base band signal. 
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A similar analysis for the quadrature channel gives two terms: 


where Im{.} is the imaginary operator and, similar to the J channel, the second term is at 
frequency 2@, and is filtered. The J and Q outputs are combined in the processor to 


form the complex signal [12]. 
C. DIGITAL BEAMFORMING 


For a simple linear array along the z axis in Figure 7, the array factor is given by: 


N N : 
AF (0,9) = ¥ w, exp| jd, cos0|= > A,e/" exp| j Bd, cos 6] (11) 
n=l 


n=1 





(n—-1)d cosé 





“*s. Wavefront 





J dy cos? 





Figure 7. _ Digital array beamforming on receive. From [12]. 


where WN is the number of elements, 8 =27/A (A is wavelength), w, = A,e/ n is a 


complex weight added by the processor and 


[2n-(N +0] , ‘i; 


n 


is the location of element n relative to the origin at the center of the array. 
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The demodulators provide spatial samples of the exponential factor in Equation 


(11): 

if =cos| jBd, cos 6] (13) 

Q,, =sin| jBd,, cos 0] (14) 
Therefore, to obtain a response equivalent to a beam scanned in the direction 0, the 
weights must be given by: 

Wy = exp|—jfd, cos 8, | (15) 


Note that compensation for any phase errors in the hardware could be included in the 


exponent of w, . They can be obtained by calibration [12]. 


D. ARRAYS WITH GROUND PLANES 


This section derives the formulas used in the Matlab codes to obtain the radiation 
pattern of a two-dimensional array with a rectangular element lattice.* Figure 8 shows a 


planar array of point sources in the x-z plane. 





Figure 8. _ Planar array centered at the origin. After [16]. 


2 This analysis is updated from unpublished notes by Prof. David C. Jenn of NPS. 
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The pattern angles 6 and ¢ and the array orientation are shown in Figure 9. 





Figure 9. _ Array orientation and pattern angles @ and @. After [17]. 


If the array is centered at the origin, then the element locations can be written as: 


2n— 1 
py Nag 


a 2 ae Tee (16) 
2m—(N, +1) 
i ag =1, »N, 


where d, and d, are the lattice dimensions and N,. and N, are the number of elements 


along x and z, respectively. The array factor is: 


N, N, 
AF = Awe sin Ocos, /BZm cos 0 (17) 


n=l m=1 


If the feeding arrangement results in a separable distribution (which is common) 


then A,,,, =A,,A, and 


Ny N, 
AF = AF. AF, => Ane Ae ante Magee cos@ (18) 


n=l m=1 
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Furthermore, if phases are introduced to scan the beam to the direction (0,,¢, ) 


then 
Ny Nz 
AF - » A,e/ Pn (sin cos ¢—sin 9, cos ¢, ) » Ae Pem (cos 8—cos 0, ) (19) 
n=1 m=1 


Finally, if the array is uniformly excited (A, 


nn =1), the sum becomes a geometric series 


with the following closed form result [18]: 


sin(N,.y,/2) sin(N,y, /2) 
sin(w,/2) — sin(w, /2) 





AF = AF, AF, = (20) 


where 


y, = fd,,(sin @ cos ¢—sin 0, cos ¢, ) 
WY, = Bd,(cos@—cos 6, ) 


Now, if the array is placed at height h above a perfect electric ground plane of 
infinite extent, then the method of images can be used to compute the field above the 
ground plane. The total array (sources plus images) can be considered as a two-element 


array along the y axis, where each element itself is a planar array as shown in Figure 10. 
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Figure 10. Two element linear array along the y axis, where each element is a planar 
array. After [16]. 


17 


The ground plane factor (GF) for a two-element linear array along the y axis, with 


the images out of phase with the sources is: 


GF = eJPhsindsing _ .— jPhsinOsing _ 9 isin ( Bhsin Osin g) (21) 
The total normalized pattern factor is obtained using the principle of pattern 
p g 


multiplication 


F(0,$) = (AF, AF) norm *GFaorm * EFnorm 
=(AF,- AF,-GF-EF) (22) 


norm 
_ sin(N,y,/2) sin(N,y,/2) 
~ N,sin(y,/2) N,sin(y, /2) 





sin (Bhsin Osin ) EFyorm 


where EF,,o;m 18 the normalized element factor. 


In general, the element factor will have both @ and ¢ components, and therefore 
F has both components as well. Referring to Figure 9, the dipoles are aligned with the x 
axis (vertical polarization). For a half-wave dipole along the x axis, the element factors 


are (see Appendix A): 


% a. 
_ Nol gr iPr cos( “sin Ocos¢) 


EF, 
2ar 1-sin? 0 cos” ¢ 


cos Ocos ¢ 





(23) 


: qT. 
itoln eo IBr cos( “sin Ocos 6) 


2ar 1-sin? Ocos* ¢ 





EFy = sing 


fifolig gr IPt 
Tr 


To normalize these equations the leading factor is simply removed. 


If the dipoles are ideal (with constant current and length, / << 1), instead of half 
wave dipoles then the terms in square brackets [ ] are reduced to 1 and / is added to the 
leading factor. Therefore, after normalization of the element factors, only the 
trigonometric functions to the right of the bracketed terms remain [14]. This concept is 


used in the Matlab codes as shown in Appendix B. 
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E. AMPLITUDE COMPARISON MONOPULSE SYSTEM 


The tracking technique for the UAV application is an amplitude comparison 
monopulse system [2]. Amplitude monopulse direction finding is a common technique 
employed in radar systems to improve the accuracy with which the direction of the arrival 
of a single pulse can be estimated. Two “squinted” beams are formed simultaneously, 


such that they overlap at the 3 dB points of the beams, as shown in Figure 11. 


-15 -10 10 15 


5 0 5 
Bearing (% 


Figure 11. Two overlapping beams. From [20]. 


By comparing the relative amplitude of the voltages of the signal being tracked 
out of the two beams, its position in the beams can be determined with an accuracy 
dependent on the signal-to-noise ratio (SNR). A common rule of thumb is that the 
accuracy of a tenth of beamwidth can be achieved with an SNR of 10 dB [18]. The 
difference of the amplitudes obtained from these two beams (the difference beam) gives 
the angular error. If the phase of the sum pattern is compared to that of the difference 


pattern, the direction of the angular error is found, as illustrated in Figure 12. 





Figure 12. Illustration of azimuth sum and difference beams. After [20]. 
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The sum beam is used for transmission and the difference and sum beams are 
used together for the receive side in monopulse systems. The sum and difference beams 


are also shown in Figure 13. 
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Figure 13. Sum and difference beam output voltages. From [20]. 


The ratio of the difference beam to the sum beam is used to generate the angle 


error difference. The tracking process uses the complex voltage ratio: 


A _ difference voltage 





(24) 
z sum voltage 
The complex voltage in Equation (25) can be separated into the real and 
imaginary terms: 
A 
Re a = Al cers (25) 
x) [z 
A 
Im . | sing (26) 
x) | 





Al. 
where 6 is the relative phase between sum and difference channels. Only the Re {3 | is 


processed, as it defines the side of the null with the sign of the ratio and the signal 


information exists only in the real term of the ratio [17]. 


20 


The plot of the pattern angle versus normalized difference to sum voltage ratio is 


shown in Figure 14. 


0.5 slope &, 
+) 
ot 
2 9 
4 
-0.5 
-1 
-15  -10 -5 10 15 


0 bs) 
Bearing (°) 


Figure 14. Normalized difference to sum voltage ratio versus the pattern angle. 
From [20]. 
The slope in the linear region creates an error signal that helps the difference null 
to correct the beam, such that the target signal is within the null. This ratio is proportional 


to the multiplication of the slope constant k, with the pattern angle in degrees as: 





6|< (27) 


where @, is the half-power beamwidth (HPBW). The principle mentioned above is 


exploited in the tracking system of this array application. It must also be noted that the 


accuracy of this techniques can be affected by unwanted multipath reflections [2]. 


F OTHER EXTERNAL CAUSES OF ERROR 


1, Multipath 


This phenomenon occurs when a signal is reflected from the external 
environment, such as the terrain and other objects. The signal reaches the receiver from 


different paths and enters the antenna from different points on its main beam or sidelobes. 


The elevation tracking errors are contributed to by smooth highly reflecting ground 
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surfaces. Diffuse scattering from buildings or rough surfaces cause both azimuth and 
elevation tracking errors. The multipath errors increase as the altitude of the target gets 


lower. 


A rms value of the multipath error is calculated from [21]: 


Oo, = 
© (8G, (peak) 


where pis the magnitude of the surface reflection coefficient, Opis antenna half-power 
beamwidth (HPBW) and G,,(peak) is the power ratio of the tracking antenna sum- 
pattern peak to the error-pattern peak sidelobe level at the angle of arrival of the image 
signal. One of the objectives in looking at a planar array is to reduce the effect of 


multipath. Adding more rows decreases the elevation HPBW (6, ) and also decreases 


O;, according to Equation (29). 


2. Crosstalk Caused by Cross-Polarized Energy 


Crosstalk, or cross-polarization coupling, is due to the cross-polarized energy of 
the source, like the undesired capacitive, inductive or conductive coupling from one 
circuit, part of a circuit or channel to another. This effect is generally negligible with a 
good array antenna element design. However, in situations where crosstalk is significant, 


a large tracking error will result in the loss of track. 
3. Troposphere Propagation 


In the troposphere, which includes altitudes between 0 through 10 km of the 
surface, the variation of the air density, dielectric constants and the thickness of the 
clouds causes beam bending. Tropospheric refraction can lead to a difference in the 
sensed direction from the bent ray and the actual direction. At standard atmosphere 


conditions and short range, this error is not significant [2]. 
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G. SUMMARY 


Normally in a digital array, each element can act as a transmitter and receiver. A 
LO is used to synchronize the elements to enable concurrent transmission and reception 
of signals. The digital beamforming concept is discussed. The ground plane effect on the 
radiation pattern due to the image created by the ground plane is derived. The tracking 
technique used is an amplitude comparison monopulse system where the sum and 
difference beams help in accurately tracking the UAV. In the next chapter, the computed 


array patterns are compared to simulated pattern using CST Microwave Studio. 
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HI. ARRAY TRADEOFF STUDY 


This chapter addresses the antenna requirements for the array design. The antenna 
requirements are the scan range, HPBW and coverage in elevation based on the UAV 


range and altitude. 
A. HALF-POWER BEAMWIDTHS 


Several array antenna analyses were done to determine the suitable array 
configuration for the tracking function. Gezer [2] designed and built a linear array of 8 
elements. The question was whether increasing the number of elements to form a planar 
array would improve the tracking and range performance. First, a 2 by 8 element array is 
compared to the original 1 by 8 element array to determine the worthiness of using a 
larger number of elements. Next, the difference and sum beam patterns are compared 
with the results generated using the CST Microwave Studio simulations from the thesis 


by Gezer. 
Simulations used the following parameters: N,=1 or 2, N,=8, d,=10 cm, 


d,=6.5 cm, h=0.25/. The azimuth and elevation patterns were computed for 2.4 GHz 


using the Matlab codes in Appendix B. Figure 15 shows the elevation patterns of arrays 


with | and 2 rows (N, =1, 2). The half-power angles are indicated by the arrows with 


the degrees shown. From Figure 15, the half-power beamwidth of a 1 row array is 


80° and the HPBW of a 2 row array is 34°. An ideal dipole is used as the array element. 
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Figure 15. Elevation plane HPBW comparison. 


In order to have the largest values of gain in the coverage region, the array can be 
tilted back such that its HPBW point touches the horizon, as shown in Figure 16. The 


range of elevation angles, as the UAV flies at an altitude of 300 m, is shown in Figure 17. 











Figure 16. ‘Tilting of array for maximum elevation coverage. 
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Figure 17. Elevation angles at various ranges. 
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The UAV is within the HPBW of a one row array for all the angles shown in 
Figure 17. For the 2 row array, the UAV is in the HPBW at ranges greater than 560 m. 
Therefore, elevation tracking would have to be done at close ranges. One advantage of 
the 2 row array is that a null is pointed at the ground. The null and the low sidelobe help 


to suppress multipath. 


The same Matlab code in Appendix B is used to generate the azimuth pattern of 
the array antenna for three different conditions, which are shown in Figure 18. The three 


conditions are as follows: 


1. lrow array at 0, =90° and ¢, =90°, 
2. 2 row array at 0, =90° and ¢, =90°, and 
3. 1 row array at 0, =40° and ¢, =90°. 


It is noticed that the azimuth pattern is the same for the first and second 
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Figure 18. Azimuth (H-plane) coverage of different arrays. 
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The elevation pattern in Figure 19 is generated for the three different 
configurations using the same Matlab codes in Appendix B. 
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Figure 19. Elevation (E-plane) coverage of different arrays. 


This tilting of the array antenna will reduce the multipath effects and the 
reflections from the earth. The null region and first sidelobe of the elevation pattern for 
the 2 row array is pointed at the ground, which reduces reflections. Due to the constraint 
on the size of the dipole cards, the closest x-direction element spacing is 10 cm. Although 
by adding 8 more elements, the signal-to-noise ratio (SNR) can be increased by a factor 
of 8, this is not an overwhelming advantage as the transmitter power of the UAV is 


sufficiently high. 


B. SUM AND DIFFERENCE BEAM MATLAB RESULTS COMPARED 
WITH CST MICROWAVE STUDIO RESULTS 


The results, using the formulas in Matlab codes for the sum and difference beams 
will be compared with those generated by CST Microwave Studio from the thesis by 
Gezer [2]. The reason for comparing the Matlab results with the CST Microwave Studio 


is to verify that the Matlab results are close to the CST Microwave Studio results, as CST 
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Microwave Studio includes the effects of the ground plane edges and mutual coupling on 
the patterns. It also has the exact dipole circuit modeled, which gives the proper element 


factor. 


CST Microwave Studio is a software package for electromagnetic analysis and 
design in the high frequency range. The Finite Integration Technique (FIT), a time- 
domain numerical method based on the integral form of Maxwell’s equations, is used in 
CST Microwave Studio. The parameters of the device that is modeled can be defined in 
detail in three dimensions. The program is especially suited to the fast efficient analysis 
and design of components like antennas, filters, transmission lines, couplers, connectors, 


printed circuit boards, resonators, etc. 


The simulation model used to run the sum and difference patterns of the array of 


dipoles with scanning in the H-plane of the array is shown in Figure 20. 





Figure 20. Top view of the array of dipoles with ground plane modeled in CST 
Microwave Studio. From [2]. 
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The Matlab codes in Appendix B are used to generate the sum and difference 
beam for comparison with the CST Microwave Studio. In order to compute the difference 


beam, AF, in Equation (18) is modified slightly. Assuming N, is even, for the 


difference beam: 


N,/2 N 
AF, = > A, ec) P2m (cos@—cosO,) _ > A, ,e/Pem (cos @—cos 0, ) (29) 
m=1 m=N,/2+1 


It can be seen that the array output of one-half of the array is subtracted from the output 


of the other half. 


The H-plane comparison results of the difference and sum beam for a scan of 0° 
from broadside are shown in Figures 21 and 22, respectively. The two methods are in 


good agreement, except at very wide angles. 
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Figure 21. Azimuth (H-plane) difference beam comparison at scan of 0° from 
broadside. 
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Figure 22. Azimuth (H-plane) sum beam comparison at scan of 0° from broadside. 


The H-plane comparison results of the difference and sum beam for a scan of 40° 


from broadside are shown in Figures 23 and 24, respectively. Again, the two methods are 


in good agreement. 
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Figure 23. Azimuth (H-plane) difference beam comparison at scan of 40° from 


broadside. 
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Figure 24. Azimuth (H-plane) sum beam comparison at scan of 40° from broadside. 


The E-plane comparison of results is shown in Figure 25. For the Matlab 
calculation, both a half-wave dipole pattern and ideal dipole pattern are shown, as 


described in Appendix B. 
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Figure 25. Elevation (E-plane) pattern comparison. 
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C. SUMMARY 


The one row array is chosen primarily because the UAV will be within the HPBW 
for all ranges. Therefore elevation tracking is not necessary. Although the SNR can be 
increased by 8 times by adding a row of 8 more elements, this will add to the cost of 
equipment and the UAV transmitter power is sufficiently high. The azimuth and 
elevation patterns of the array were examined. The sum and difference beams for the 
Matlab simulations are in close agreement with the CST Microwave Studio simulations. 
The printed circuit dipole element was found to have an E-plane pattern HPBW between 


that of an ideal dipole and a half-wave dipole. 
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IV. FREQUENCY MODULATION (FM) DEMODULATOR 


A. INTRODUCTION 


Apart from tracking the UAV, there is a requirement to recover the information 
from the signals transmitted by the UAV. The signal of interest will be either a NTSC 
video signal or audio signal, that is frequency modulated (FM) onto a 2.4 GHz carrier. 


The modulation characteristics of the FM video and audio signals are shown in Table 1. 


Parameters | Bandwidth | FM Modulation 
Modulation | Index 


Video and | 6 MHz Wideband >1.0 
Audio 


Table 1. Parameters of FM video and audio signals. From [22]. 





UAV surveillance has attempted to take advantage of cheap wireless video 
transmitters or receivers that operate at 2.4 GHz in the ISM band (license free band). 
Typical low power transmitters have power ratings in the range of 10-100 mW, although 
more powerful devices are available. For UAV applications, longer ranges are possible 
using higher power. The 2.4 GHz wireless video transmitters, which can transmit a 
broadcast quality picture, generally use a FM modulation method for audio/video 
transmission. One thing to note is that 2.4 GHz wireless video systems are not the only 
equipment that uses this frequency band. For example WLAN systems and microwave 
ovens operate at the same frequency band and can cause interference in the operation of a 


2.4 GHz video link. 


The output of each array element is an FM signal with a 2.4 GHz carrier. These 
signals need to be weighted and summed for beamforming, and then demodulated to get 
the video or audio. The beamforming can be done at the carrier frequency or at baseband. 


Either way, it is necessary that the incident wave phase change between elements is 
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retained for beamforming. There are several approaches to demodulating the FM signal at 


the antenna. Table 2 shows the types of FM detectors that are used for the demodulation 


of FM signals. 


Name of Type of Circuit used 
detector detector 


Frequency Non-coherent | Quadrature demodulator (Bilotti) 
discriminator Foster-Seeley detector 
Balanced frequency discriminator 


Phase-locked | Coherent Analog phased-locked loop 
loop (PLL) 

Digital Coherent Analog-to-Digital Converter 
sampling (ADC) 


Table 2. Types of FM detectors. From [24] 





The ADC for demodulation of the FM signals (also known as the digital 
demodulator), at some point, converts the analog signal into a digital representation prior 
to demodulation. The most direct approach is to sample the received signal out of the 
antenna, but this requires a very high-speed ADC. If ten times oversampling is used 
(which is common), the sampling frequency would have to be 10x 2.4 GHz = 24 GHz. 
The fastest ADCs today are in the 2-4 GHz range. Two common approaches to 
implement this are to sample at a lower intermediate frequency (IF), e.g., around 30 MHz, 
or to mix the signal down to baseband, where I/Q components are required to retain all 
the information in the signal. The Foster-Seely detector and balanced frequency 
discriminator convert FM into amplitude modulation (AM) and use an AM detector to 
recover the message signal from the converted waveform. The quadrature demodulator 
converts FM into both phase modulation (PM) and AM and either type of detector can be 


used to recover the message signal. This is referred to as Bilotti’s method. 


A local oscillator (LO) signal is necessary for downconverting from 2.4 GHz to 


baseband. The LO signal can be provided by an oscillator at the antenna, or by extracting 
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(recovering) the carrier from the received signal. Each of these methods has advantages 
and disadvantages. A phase-locked loop (PLL) is one way to recover the carrier. A PLL 
is an electronic control system that generates a signal that has a fixed relation to the phase 
of a "reference" signal. A phase-locked loop circuit responds to both the frequency and 
the phase of the input signals, automatically raising or lowering the frequency of a 
controlled oscillator until it is matched to the reference in both frequency and phase [25]. 
If the analog signal processing, prior to sampling into the digital domain, is done with 
phase-locked sources, the digital demodulator can support either coherent or non- 


coherent demodulation. This is a key potential benefit of the digital solution [26]. 


Several of the methods for the demodulation of the FM signals will be examined 


to assess their applicability to the tracking array. 
B. LOCAL OSCILLATOR AND DEMODULATION APPROACHES 


Three FM demodulation methods are briefly explained in this section. A closely 
related issue is how the local oscillator signal is generated. The LO signal can be 
extracted from the incoming signal or generated locally. Both of these methods may 


require feedback to keep the LO locked on the incoming carrier. 
1. Phase Locked Loop 


A PLL is a control circuit that is capable of acquiring a signal and locking on to it. 
The PLL system consists of three parts: phase comparator, low-pass filter (LPF), and 
voltage-controlled oscillator (VCO). All parts are connected to form a closed-loop 


frequency-feedback system, as shown in Figure 26. 
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Figure 26. PLL Block Diagram. From [27]. 
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With no signal input applied to the PLL system, the error voltage at the output of 


the phase comparator is zero. The voltage, V(t) , from the LPF also is zero, which causes 
the VCO to operate at a set frequency, f,, called the center frequency. When an input 


signal is applied to the PLL, the phase comparator compares the phase and frequency of 
the signal input with the VCO frequency. It then generates an error voltage proportional 
to the phase and frequency difference of the input signal and the VCO. The error voltage, 


V(t), 1s filtered and applied to the control input of the VCO. V(t) varies in a direction 


that reduces the frequency difference between the VCO and signal-input frequency. 
When the input frequency is sufficiently close to the VCO frequency, the closed-loop 
nature of the PLL forces the VCO to lock in frequency with the signal input; 1.e., when 
the PLL is in lock, the VCO frequency is identical to the signal input, except for a finite 
phase difference. The range of frequencies over which the PLL can maintain this locked 
condition is defined as the lock range of the system. The lock range always is larger than 
the band of frequencies over which the PLL can acquire a locked condition with the 
signal input. This latter band of frequencies is defined as the capture range of the PLL 


system [27]. 
2. Local Oscillator With No Feedback 


Using an independent LO with no feedback is similar to a PLL in the free running 
mode, as shown in Figure 27. Ideally the VCO frequency should be exactly the same as 


the carrier frequency, f.. If not, then the reference for frequency demodulation will be 


off. If the carrier frequency of the transmitter is stable and accurate and the VCO 
frequency is also stable and accurate, then the baseband signal can be demodulated. In the 
next two sections, Matlab simulations are demonstrated to estimate how much frequency 


error between the input signal and VCO can be tolerated. 


Phase 
Comparator 


Figure 27. Stand-alone VCO for a local oscillator. 
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Signal 


Input 
Signal 





3. Bilotti’s Quadrature Demodulator? 


The modulation of an FM signal is the deviation of a carrier from its nominal 
frequency. To demodulate this signal, its frequency deviation is converted to either phase 
or amplitude and the change in phase or amplitude is detected. In the quadrature 
demodulator (also called Bilotti’s method), the modulated carrier is passed through an LC 
tank circuit that shifts the signal by 90° at the center frequency. The direction of 
frequency deviation will cause the phase shift to be either greater or smaller than 90°. 
The LC tank circuit simultaneously converts the FM to amplitude modulation (AM), so 


an amplitude detector could also be used. 


Bilotti's quadrature demodulator uses a phase shift network and a phase detector, 


as shown in Figure 28 [29]. 





Figure 28. Quadrature demodulator block diagram. From [29]. 


The signal out of the antenna is v;,. The phase detector compares the phase of the 
input signal (1, ) to vz, the signal generated by passing v, through a phase shift network. 
This phase shift network includes an LC tank (L,R, and C p) and a series reactance 
(C, ). The network gives a frequency-sensitive phase shift. The phase detector discussed 


is the bipolar, double-balanced multiplier popularized by Bilotti. The output of the 
multiplier (7, ) is filtered, which results in a DC level that changes as the input frequency 


changes. 


3 Much of this material in section A and B is from [29]. 
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To derive the transfer function of the quadrature demodulator, the phase shift 


network is drawn as a small-signal circuit model in Figure 29. 





Figure 29. Small signal model of the quadrature phase-shift network. From [29]. 


The impedance (Z,, ) of the parallel combination of L, R, and C,, is: 


ROL 
Z,(s)= 2 (30) 


= 
CpRpLs +Ls+Ry 
where s= jo. The ratio of vy over y is the ratio of impedances Z,(s) over 


(Z,(s)+1/sC,). Simplifying this ratio gives: 


V2 _ Cis 





- i G1) 
MW (OSHC Nis tail 
P s R 
Dp 


Equation (31) is used in the Matlab code in Appendix C to generate the plot of the phase 


angle versus the frequency. The resonant frequency @, of this filter is: 


On = 2 (32) 


" {L(C,+Cp) 
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Cc. SIMULATION RESULTS 


In this section, simulation results are presented for several FM demodulation 


techniques. 
1. Quadrature Detector (Bilotti’s Method) 


In Bilotti’s quadrature detector, the output can be applied to either an AM or PM 
detector. In the present application, the mixer and filter in Figure 28 are incorporated in 
the AD8347. The output i, is in I/Q format. The demodulation is done in a computer 
using the sampled I/Q data. Therefore, either AM or PM could be used to recover the 
video signal. 

Equation (31) and the Digikey catalog [30] are used to compute and select the 


appropriate values for the capacitances C p and C, and inductance L to obtain a carrier 
frequency of 2.4 GHz. Using off the shelf values of C P and C,, both are 5 pF. From 


available off the shelf values of L, 45 nH gives the exact carrier frequency of 2372.542 


MHz. These lumped element values are used in the simulations. 


Using the Matlab codes in Appendix C and varying the values of R,, various 


plots for the 6 MHz video bandwidth are obtained as shown from Figures 30 to 33. 
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Figure 33. 


It can be observed that as R 


p increases, more phase angles can be covered within 


the 6 MHz bandwidth about the carrier frequency of 2372 MHz. However, the plots 


become more and more non-linear. The characteristics are summarized in Table 3 for the 


6 MHz bandwidth at about 2372 MHz. 


Angle range,° 


Behaviour 





80 to 104 


Linear 





71 to 114 


Linear 





64 to 123 


Linear 





51 to 140 


Non-linear 





40 to 150 


Non-linear 





29 to 157 


Non-linear 





20 to 166 


Non-linear 











Table 3. | Summary of the plots in Figures 30 to 33. 


If the computer processor uses AM, then a large dynamic range in amplitude is 
best. However, if PM is used then a linear phase curve over a large range of angles is best. 


various 


Using the Matlab codes in Appendix C and varying the values of R,, 


plots for the 22 kHz audio bandwidth are obtained, as shown in Figures 34 to 37. 
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It is noted that as Rp 


cover a larger range of phases, but the non-linear characteristic becomes more obvious. 


For the 22 kHz bandwidth, R P is chosen to be 1 MQ. 


The quadrature circuit could be implemented in the digital array, as shown in 
Figure 38. The reference signal can be obtained using a special auxiliary antenna, or it 


could be tapped off of one of the elements. However it is done, the same LO source must 


be used for all elements. 
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increases from 510 kQ to 4.04 MQ, the higher values 
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Digital array using a quadrature detector. 


Figure 39 shows a Matlab Simulink model for the array in Figure 38, up to the 


computer block. A two-tone signal, one sine wave at 900 Hz and another at 1700 Hz, is 


used as a source that is converted to FM. The simulation operates at 50 kHz, rather than 


2.4 GHz, to greatly reduce the simulation time. The frequency deviation of the FM is set 


at 1000 Hz, which gives a 0.5 modulation index, assuming a maximum modulating 


frequency of an FM signal is 2000 Hz. A sampling frequency, which is 10 times the 


simulation frequency, is used to obtain a more refined simulation result. The same 


sampling frequency will be used for subsequent models in this section. The RLC circuit 


is modeled by its transfer function, and its resonant frequency was set to 50 kHz. 
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Scope 


Simulink model of Bilotti’s quadrature demodulator. 


The FFT of the source and recovered signals are shown in Figure 40. The 
frequency components of the two-tone signal are clearly visible in the recovered signal. 


This shows that the source signal can be recovered using Bilotti’s method. 
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Figure 40. FFT of the recovered and source signals in Bilotti’s simulink model. 


2. Local Oscillator With No Feedback 


Figure 41 shows the block diagram for an array with an independent LO. The LO 
will synchronize the AD8347 demodulators. As this is a method without feedback, 
simulations are done to consider the effects of a phase drift in the LO. The raw I and Q 
components from the AD8347 demodulators are processed to perform the beamforming 
and demodulate the FM video signal. The carrier frequency of 50 kHz with a 1 kHz 
deviation is used in the simulation for the modulator for transmitting the two-tone signals. 


Figure 42 is the Simulink model of Figure 41. 
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Figure 41. Open loop digital quadrature detector. 
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Figure 42. Simulink model of open loop digital quadrature detector. 


Figure 43 shows the input and demodulated signals in the time domain. Figure 44 
shows their spectra in the frequency domain. The recovered signal matches the source 


signal in both the time scope and the FFT spectrum scope. 
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Figure 43. Time scope of the open loop digital quadrature detector in Simulink. 






Source signal 
in FFT spectrum 


Magnitude-squared, dB 


-100 


0 2 25 3 35 4 45 5 
20 
10 Recovered signal 


in FFT spectrum 


Magnitude-squared, dB 
8 








2. 35 4 45 5 


2 5 3 
Frequency, kHz 


0 05 1 1.5 


Figure 44. FFT spectra of the source and recovered signals. 
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3. Digital PLL Quadrature Detector 


Figure 45 shows how the PLL is integrated into the digital array to perform the 
beamforming and the FM demodulation. The output I and Q components are converted to 
the digital format to be processed in the computer. In the computer, the digital I and Q 
signals are low pass filtered before performing the beamforming and FM demodulation 
The error signal is computed and fed to the variable control oscillator (VCO) to adjust the 
VCO frequency. The same PLL is used for mixing all the elements from the other 
channels. 
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Figure 45. Digital PLL quadrature detector. 


The Simulink model in Figure 46 was used to simulate the digital PLL quadrature 
detector in Figure 45. The carrier frequency of 50 kHz with a | kHz deviation is used for 
the modulator for transmitting the two-tone signals. The VCO with the following 
parameters is inserted in the feedback loop: output amplitude of 3 V, quiescent frequency 


of 50 kHz, input sensitivity of 500 Hz/V and a sampling frequency of 500 kHz. 
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Figure 46. Simulink model of digital PLL quadrature detector. 


Figure 47 shows the input and demodulated signals in the time domain. Figure 48 
shows their spectra in the frequency domain. It is noticed that the source signal matches 
the recovered signal for the digital PLL quadrature detector. 
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Figure 47. Time scope of the digital PLL quadrature detector in Simulink. 
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Figure 48. FFT spectrum of the digital PLL quadrature detector in Simulink. 


D. PHASE ERRORS AND LO FREQUENCY DRIFT 


There is a need to study the effects of the frequency drift in the LO and the 
relative phase errors in the array channels. The relative phase errors cause beamforming 
errors, which lead to tracking errors. In this section, the effects of phase errors of a 4- 
channel array modeled in Simulink are shown. The frequency drift in the LO is also 


simulated. 
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1. 


Relative Phase Errors in the Array Channels 


The Simulink block diagram designed to model the phase errors in the array 


channels is shown in Figure 49 below. The carrier frequency of 50 kHz with a | kHz 


deviation is used in the simulation for the modulator for transmitting the two-tone signals. 


Phase shifters are introduced to the open loop digital quadrature detector Simulink model. 


Different delay line sizes are introduced to the other three channels to model the effects 


of phase shifters on the incoming signal. Ideally, the delays cause the sum output to have 


a maximum at the signal angle of incidence. The Simulink results, shown in Figure 50, 


indicate that the signal could be recovered when phase shifter errors are introduced. 
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Simulink model for phase errors simulation analysis. 
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Figure 50. Simulation results for phase errors in a 4-channel array. 


2. LO Frequency Drift 


The frequency drift in the LO is simulated in the Simulink model, shown in 
Figure 51, at a carrier frequency of 5 MHz to analyze the effects to the recovered signal. 


A frequency offset of 100 Hz is introduced to the transmitter frequency. 
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Figure 51. Simulink model for LO frequency drift analysis. 


The simulation result is shown in Figure 52. 
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Figure 52. Simulation result at 5 MHz with 100 Hz frequency drift. 
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E. SUMMARY 


Simulation tools are developed for array analysis of the three FM demodulation 
techniques: Bilotti’s method, local LO with no feedback and digital phase locked loop 


quadrature detector. 


The LO with no feedback is preferred as it uses the least circuitry with 
satisfactory performance for the recovery of the two-tone signal. Beamforming can also 
be performed easily, as the raw I and Q data is preserved. This technique also eases the 
computer processing. Relative phase errors in the array channels and LO frequency drift 
are introduced into the Simulink model. The simulation results indicate that the recovery 
of the two-tone signal can still be achieved. Hence, the LO with no feedback technique 


can be used. 
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Vv. SYSTEM ARCHITECTURE 


This chapter discusses the hardware necessary to build the tracking array. To keep 


the cost down, commercial-off-the-shelf hardware is used. 


DESCRIPTION OF TRACKING ARRAY HARDWARE 


A general block diagram of a digital receiver channel is shown in Figure 53. The 


received FM video or audio signal undergoes the following processes: 


1. Amplified by the low noise amplifier (LNA) and then fed to the RF input 
of an AD8347 quadrature demodulator. 


2. Downconverted to baseband by the AD8347 quadrature demodulator, 
operating from 800 to 2700 MHz with a modulation bandwidth of about 50 
MHz. The output signals are in I/Q format. 


3. Sampled, filtered and then sent to the processor. 


4. Processing of the detected signals is done to form the beam by combining 
the I/Q signals from all receive channels. The FM is removed and the video 


signal is recovered. 


5. A common local oscillator is used to synchronize the devices to avoid 


phase and drifting errors between the elements. 


These processes are illustrated in Figure 53. 
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Figure 53. Array architecture for one channel. 
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The dipoles, LNA, demodulator, local oscillator, power supplies, power 
adjustment circuit and signal processing are discussed briefly, and follow the description 


in Gezer [2]. 
iE Dipoles 


The dipoles were designed in the CST Microwave Studio by Professor David C. 
Jenn and fabricated by Cirexx Corporation. They operate at 2.4 GHz. The return loss of 
the dipoles was measured and found to be at least 15 dB. Reference [8] discusses the 


dipole element return losses more in detail. 
2 Low Noise Amplifier (LNA) 


The LNA amplifies the weak signals captured by the dipoles. The LNA has a low 
noise figure (2.0 dB) and high gain (24 dB), which improves the performance of the array. 
Because the overall noise figure is dominated by the first few stages, the total noise of the 
cascaded stages is reduced by the LNA. Hence, the fact that the LNA has a very low 


noise figure is a key factor in signal detection. 


The LNA must be located very close to the dipole in order to minimize the losses 
in the feed line. For our setup, the LNAs are connected to the dipoles with a barrel to 
minimize the cable loss. This SNR will be verified later in Section B. The LNA-2700 by 
RF Bay Inc. (shown in Figure 54) was selected as the LNA because of its reasonable 


price as well as its low noise figure. 





Figure 54. LNA 2700 by RF Bay Inc. From [2]. 
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3. Demodulator 


The demodulator downconverts the signal from the carrier frequency to the 
baseband. The demodulator selected for the array is the AD8347 Quadrature 
Demodulator manufactured by Analog Devices. Figure 5 shows one of these 
demodulators used in the system. This COTS demodulator operates within the frequency 
range 0.8 through 2.7 GHz and requires supply voltages between 2.7 and 5.5 V. The 
voltage VGIN determines the amplification of the received signal. Reference [31] 


includes a very detailed spec sheet of the AD8347 EVAL Quadrature Demodulator. 


4. Local Oscillator 


The LO is designed and manufactured by Z-Communications. The LO signal at 
2.4 GHz is generated by a V800MEI10 oscillator on the board and fed into the 


demodulators. 


5. Power Supplies and Power Adjustment Circuit 


The two linear power supplies used for the array are model HP 6255 A by 
Hewlett Packard. The power supplies provide +5 and +12 V. The +12 V is applied to 
VGIN of the demodulator after being adjusted with a VGIN voltage regulator, designed 
by Bob Broadston from NPS. The LO and LNAs are provided with +12 V by the power 
supply, which has maximum currents of 6.00 A, 2.50 A for the voltages of 5 V and 12 V, 


respectively. 


6. Analog-to-digital Converters 


The ADCs sample the baseband I/Q signals out of the AD8347. The minimum 
sampling rate is set by the baseband bandwidth, which is approximately 6 MHz for the 
FM video signal. For 5 times oversampling, this would be 60 MHz. For the audio signal, 
it would be about 110 kHz. 


Candidate hardware for the audio case is the National Instruments Compact 


Reconfigurable I/O (cRIO) product line. The reconfigurable chassis is the heart of the NI 
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cRIO system and contains a Field Programmable Gate Array (FPGA). The cRIO core is 
powered by a user-programmable FPGA configured by a special software called 
LabVIEW, that serves as a high-speed interface to the I/O modules. The LabVIEW 
FPGA Module is used to program the FPGA on the cRIO reconfigurable chassis. FPGAs 
are silicon chips with unconnected logic gates. The flexible design of FPGA enables the 
user to target any device compatible with NI cRIO and to customize hardware functions 
for its system. It is possible to execute multiple tasks simultaneously using the LabVIEW 


FPGA module. 


NI provides different I/O modules for different purposes. The four analog output 
signals from the demodulator are input to an analog input (A/I) module to be sampled. 
They are then converted into digital signals for signal processing and analysis. The A/I 
module selected for the audio case is the NI cRIO-9215. The NI cRIO-9215, shown in 
Figure 55, provides 4 simultaneously sampled analog inputs with a rate of 100 kS and a 


resolution of 16 bits. 





Figure 55. NI cRIO-9215 Analog Input Module. From [32]. 


The NI cRIO-9215 analog input modules are connected directly to the FPGA 
hardware. This FPGA hardware is the reconfigurable chassis NI cRIO-9104 with the real- 
time embedded controller cRIO-9004. The cRIO-9104 chassis has 8 slots and a RIO 
FPGA core with a RAM of 196 kB. The FPGA core has an individual connection to each 
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NI cRIO-9215 analog input module and can read information from each module. The 


module operation is precisely synchronized with 25 ns resolution. 
B. RECEIVER MODULE ANALYSIS 


This section presents an estimate of the operating range based on the minimum 


power required by the demodulator (AD8347). 
1. Sensitivity and Dynamic Range 


The maximum RF power that is allowed at the receiver input (without saturating 
it) is limited by the demodulator’s 1 dB compression point (IP1dB). From [9], the 
demodulator’s IPIdB ~ -—22 dBm with no Automatic Gain Control (AGC). With 
reference to Figure 53, the LNA gain is 24 dB and the cable losses on the two sides of the 
LNA are 0.2 dB and 0.4 dB. The maximum RF power that is allowable to the receiver 


antenna is: 


r. 


Po =-22+0.2 -24+0.4= —45.4 dBm 
P. gives the closest range without over driving the demodulator. 


The minimum RF power that is required at the receiver module for operation 
depends on the demodulator’s sensitivity. From [9], the demodulator mixer output 
voltage is linear with RF input power from —65 dBm to —5 dBm. Hence, with reference 


to Figure 53, the minimum RF power at the receiver antenna is: 


P 


Tin 


—65 + 0.2 —24 +0.4 = —88.4 dBm 


—118.4 dBW = 1.445x 10’? W (33) 


Hence, the dynamic range of the receiver module is —45.4— (— 88.4) = 43 dB. 
2. Maximum Range 


Consider the distributed array block diagram in Figure 56 where 
R, =the distance between the UAV and the n" receive element, 
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P =transmitted power from the UAV, 


th 


P. =received power at n” element, 


G, = gain of transmitter antenna of UAV, and 
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Module oe] 
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N receive elements 2 tt 





Figure 56. Distributed array of receiver antennas. 


At any one time, N elements are receiving signals from the UAV, where WN is the 
total number of elements in the array, which in our case is N = 8. From Figure 56, the 
one-way Friis equation is used to determine the maximum operating range, R,,,.. It is 
considered that there is no polarization mismatch, input mismatch or other system losses. 


The one-way signal strength from a transmitter to a receiver element is: 


2 
FE = BEA (34) 
(47) R,” 


n 


To include receiver sensitivity in Equation (34), P, is used to determine the 


maximum range: 


(35) 
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It will be assumed that P 


- is 100 mW in this case. The gain of each of the receiver 


elements is G, and calculated to be 6.3 = 8 dBi. Thus, substitute P = 100 mW, G, = 1.6, 


G, =6.3, and P. = —88.4 dBm into Equation (35) to get R,,, = 8.31 km. 


C; SUMMARY 


A description of tracking array hardware was the subject of this chapter. The 
possible hardware components were briefly introduced. From the block diagram of the 
receiver module, the minimum power for the receiver module was used to determine the 


maximum operating range of the array using the one-way Friis equation. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY AND CONCLUSIONS 


As a part of an ongoing project, this research’s aim was to design a digital phased 
array for receiving and tracking a 2.4 GHz FM video signal from an unmanned air 
vehicle. Tracking the signal from the UAV will establish a data link for data transfer such 


as video, image or audio between the UAV and the ground station. 


An atray tradeoff study was done to determine the benefits of the 1 row and 2 row 
arrays with 8 elements per row in terms of HPBWSs of the azimuth and elevation patterns. 
The one row array was chosen primarily due to its wide elevation pattern coverage. This 
eliminates the need for elevation tracking, as the UAV is within the HPBW for all UAV 
ranges of interest. Adding a row of 8 more elements to increase the SNR will increase the 
cost of equipment. The additional SNR is not required, as the UAV transmitter power is 
sufficiently high. The sum and difference beams for the Matlab simulation, based on the 
derived pattern function of the one row array, are in agreement with the CST Microwave 
Studio simulations from the previous project by Gezer [2]. The printed circuit dipole 
element was found to have an E-plane pattern HPBW (80°) between that of an ideal 
dipole (90°) and a half-wave dipole (78° ). 


Several approaches for the demodulation of FM signals were explored. The 
considerations for the choice of the approaches for the FM demodulation are minimum 
circuitry, the ability to perform beamforming, recovery of the transmitted signal and the 
susceptibility to phase errors and LO drift. Three FM demodulation techniques are 
considered: Bilotti’s quadrature demodulator, a local oscillator with no feedback and a 
digital phase locked loop quadrature detector. All three use digital beamforming and 


signal processing to demodulate the FM. 


The Simulink models of the three FM demodulation techniques show satisfactory 
performance for the recovery of a two-tone signal. The simulations are done at a carrier 


frequency of 50 kHz, instead of 2.4 GHz, to greatly save on the simulation time. A LO 
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with no feedback is chosen, as it uses the least circuitry with satisfactory performance for 
the recovery of the two-tone signal. Beamforming can also be performed easily, as the 
raw I and Q data are preserved. This technique also eases the computer processing. 
Simulations were done with phase errors introduced in the channels and a frequency 
offset in the LO without feedback. The Simulink results indicate that the recovery of the 


two-tone signal can still be achieved with the expected LO drift. 


Finally, the possible hardware components are briefly introduced in Chapter V. 
From the block diagram of the receiver module, the minimum power for the receiver 
module was used to determine the maximum operating range of the array using the one- 


way Friis equation. 
B. RECOMMENDATIONS FOR FUTURE WORK 


Simulations and studies should be conducted to analyze the following effects on 


the array and demodulation of FM signals. 
1. Implementation of Signal Processing 


The software necessary to do the beamforming and FM demodulation must be 
acquired or developed. The beamforming and tracking are relatively straightforward, and 
have been demonstrated by Gezer [2]. However, an acquisition method must be devised 
to acquire the signal when the antenna and processor are first turned on (or when the 
transmitter signal first appears). If the PLL method is used, the processor must compute 


the VCO error signal. 
2. Hardware Tradeoff Studies 


An appropriate ADC must be selected and interfaced with the AD8347 
demodulator board, and controller software must be written. The hardware must be 
selected to minimize the amount of processing and control commands, so that real-time 


operation is possible. 
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APPENDIX A. DIPOLES OF ARBITRARY ORIENTATION! 


The dipole is centered at the origin as shown in Figure Al. The dipole direction 


is given by the vector: 


[=xsin0, cos¢, + ysinG, sing, + Zcos 0, 





u Vv 


a a Ww, 


a 
The dipole direction cosines are: 
x direction cosine: u, = cosa, =sin@, cos ¢, 


y direction cosine: vy = cos @,, = sin 6, sin g, 


z direction cosine: w, =Ccos@, = cos@, 


The length of the dipole is L, and for a half-wave dipole the current distribution is given 
by the following function of the path length variable @’ : 
I(l)=1, cos(ké’), -L/2<0'<L/2 


where J, is the maximum value at the feed point. For an ideal (Hertzian) dipole with 


constant current: 
MO.= Ey AALS S112, Len. 
The direction cosines for the far field observation point are: 


x direction cosine: u = sin @cos @ 
y direction cosine: v = sin @sin ¢ 
z direction cosine: w=cos@ 


4 The materials for Appendix A is from “Arbitary dipoles” unpublished notes of David C. Jenn. 
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Figure Al: Dipole of arbitrary orientation centered at the origin 


The far-field radiation integral can be written as: 


ane ry 65 SCL/2 ee ae 
E(6,@) = LE to! eit f I(0')eIBE (ef) ap) 
ae 5/2 


where it is understood that the 7 component is discarded. For the ideal dipole the 


integral is: 


L/2 ieee 
i 1eyeibelreD gy = InLsin| 
-L/2 


BLiet 





}. Inl (L<<A) 


and the resulting spherical field components are: 
Eg — £0, 4) 06 = IP MolmE 6 iPr? 66 
Arr 


Es ~ £0, 4)0g = —IBtolm” .- i8r 7, g 
4ar 
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From the transform tables: 


6 = &cos Ocos¢ + $cos Osing — 2 sin O 
g =—ksind + fcos¢ 


For a Hertzian dipole along the z axis: (0, =¢, =0°,f =2) 


and the fields reduce to the familiar results: 


_ iBitolnk.-ipr 
4ar 


Eo sin 0 
Ey =0 


For a half-wave dipole (L = 4/2) the integral is: 





iu L/2 
ie I, cos( Ble jBt( PD) apr Leto (Fei ul [Bl (Fef)cos( Bl’) + Bsin( Be’) | 
-L/2 [pe(ref)] +p? , 


After applying some trigonometric identities, the final result for the far electric field is: 


. s —jBr cos(ZFe/] 
cc 2ar a\2 
He 1—-(fef) 








2ar sind 


_ ily "onl see 


Ey =0 
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For a half-wave dipole along the x axis: 0, =90°,g, =0°,f=£ 


Pel =sinOcos¢=u 
def =cosOcos¢ 
gef =-sing 


‘ WT .. 
iol n Pet cos“ sin dcos¢) 


ul 5 x [eos Acos¢@ 
2ar 1—sin“ Ocos* ¢ 


Eg = 


; —jfr cos(% sin Acos 6) 
E,=—wteime 2 sind 


2ar 1-sin? 0 cos” ¢ 
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APPENDIX B. MATLAB CODE FOR AN ARRAY OVER A 
GROUND PLANE 


This Matlab code was used to compute the array pattern for | row and 2 row arrays and the sum 
and difference beams, using the formulas derived in Chapter II. 


linear array of dipoles with a ground plane 
% 7 is array axis; ground plane is the xz plane 
% y is normal to ground plane 
% option to select dipoles along x or z 
% fixed to include Etheta and Ephi 
% two dimensional array (scans in both x and z dimension) 
% half-wave dipoles 
clear 
warning off 
rad=pi/180; 
f=2.4e9; 
wave=3e8/f; 
beta=2*pi/wave; 
dz=0.065; 
Nz=8; %no. of elements in a row 
dx=0.1; 
Nx=1; %To define the no. of rows 
h=0.25*wave; 
% select the dipole direction: 'x' or 'z' 
dipdir='x'; % dipoles x directed (parallel case) 
%dipdir='z'; % dipoles z directed (collinear case) 
% select the dipole type: 'h'=halfwave, 's'=short (constant current) 
diptype='s'; 
%diptype="'h'; 
%phi=0 cut 
% phi=0 vary theta (x-z plane) 
phi=90*rad; 
theta=[0:1:180]*rad; w=cos(theta); 
thetas=90; phis=;%theta defines the yz plane; phi defines the xy plane; 
% x is along the direction of the short dipole; 
%y is perpendicular to the ground; 
%z is perpendicular to the short dipole and parallel to the ground plane. 
cts=cos(thetas*rad); u=sin(theta)*cos(phi); v=sin(theta)*sin(phi); 
cps=cos(phis*rad); 
alphaz=-beta*dz*cts; 
alphax=-beta*dx*cps; 
psiz=beta*dz*w+alphaz; 
psix=beta*dx*u+alphax; 
if dipdir=='z' 9% colinear dipoles (z directed) 
if diptype=='s' 
EFtheta=2*sqrt(1-w.*2).*sin(beta*h*v); 
EFphi=0; 
end 
if diftype=='h' 
EFtheta=2*sqrt(1-w.*2).*sin(beta*h*v).*cos(pi*w/2)./sqrt(1-w.*2+1e-5); 


de 


EFphi=0; 
end 
end 
if dipdir=="x' % parallel case (x directed) 
if diptype=='s' 
EFtheta=2*cos(phi).*w.*sin(beta*h*v); 
EFphi=-2*sin(phi).*sin(beta*h*v); 
end 
if diptype=='h' 
EFtheta=2*cos(phi).*w.*sin(beta*h*v).*cos(pi*u/2)./(1-u.42+ Le-5); 
EFphi=-2*sin(phi).*sin(beta*h*v).*cos(pi*u/2); 
end 
end 
AF=sin(Nx*psix/2)./sin(psix/2).*sin(Nz*psiz/2)./sin(psiz/2); 
caftheta=abs(AF.*EFtheta)/2/Nx/Nz; 
cafphi=abs(AF.*EFphi)/2/Nx/Nz; 
cafthdb=20*log 10(caftheta+ le-5); 
cafphdb=20* log 10(cafphi+1e-5); 
caf=sqrt(cafphi.’2+caftheta.”2); 
cafdb=10*log10(caf.’2); 
figure(1),clf 
subplot(121) 
polardb(theta/rad,cafthdb, '-') 
title({'Polar dB plot E_\theta, \phi=',;num2str(phi/rad),'*o']) 
subplot(122) 
polardb(theta/rad,cafphdb,'-') 
title(['Polar dB plot E_\phi, \phi=',num2str(phi/rad),'“o']) 
hold off 
%theta=90 cut 
clear theta phi AF caf cafthdb cafphdb EFth EFph psi u v w cts ct 
% theta=90 vary theta (x-y plane) 
phi=[0:1:180]*rad; 
theta=89.99*rad; ct=cos(theta); 
% note to always mantain the cut theta at 89.99 or 90.01 
u=sin(theta)*cos(phi); 
cts=cos(thetas*rad); v=sin(theta).*sin(phi); w=ones(1,length(v))*ct; 
alphaz=-beta*dz*cts; 
psiz=(beta*dz*cos(theta)+alphaz)*ones(length(phi), 1)’; 
psix=(beta*dx*u+alphax); 
if dipdir=='z' % colinear dipoles (z directed) 
if diptype=='s' 
EFth=2*sqrt(1-w.2).*sin(beta*h*v); 
EFph=0; 
end 
if diftype=='h' 
EFth=2*sqrt(1-w.42).*sin(beta*h*v).*cos(pi*w/2)./sqrt(1-w.*2+1e-5); 
EFph=0; 
end 
end 
if dipdir=="'x' % parallel case (x directed) 
if diptype=='s' 
EFth=2*cos(phi).*w.*sin(beta*h*v); 
EFph=-2*sin(phi).*sin(beta*h*v); 
end 
if diptype=='h' 
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EFth=2*cos(phi).*w.*sin(beta*h*v).*cos(pi*u/2)./(1-u.42+1e-5); 
EFph=-2*sin(phi).*sin(beta*h*v).*cos(pi*u/2); 
end 
end 
AF=sin(Nx*psix/2)./sin(psix/2).*sin(Nz*psiz/2)./sin(psiz/2); 
caftheta=abs(AF.*EFth)/2/Nx/Nz; 
cafphi=abs(AF.*EFph)/2/Nx/Nz; 
cafthdb=20*log10(caftheta+le-5); 
cafphdb=20* log 10(cafphi+1e-5); 
caf=sqrt(cafphi.’2+caftheta.’2); 
cafdb=10*log10(caf.’2); 
figure(2),clf 
subplot(121) 
polardb(phi/rad,cafthdb, -') 
title({'Polar dB plot E_\theta, \theta=';num2str(theta/rad),'“o']) 
subplot(122) 
polardb(phi/rad,cafphdb, '-')% plot for E phi figure 2 
title({'Polar dB plot E_\phi, \theta=',num2str(theta/rad),’“o']) 
hold off 
% clear theta phi AF caf cafdb EFth EFph psi 
% 3-D plot % 
iplt3d=1; % skip 3d plots if =1 
if iplt3d==0 
it=0; 
for th=0:3:180 
jt=0; it=it+1; 
for ph=0:3:360 
jt=jt+1; 
thr=th*rad; 
phr=ph*rad; 
phd (it,jt)=ph; 
thd(it,jt)=th; 
ct=cos(thr); 
st=sin(thr); 
u=sin(thr)*cos(phr); 
v=sin(thr)*sin(phr); 
w=cos(thr); 
xsiz=beta*dz*cos(thr)-alpha; 
xsix=beta*dx*u; 
E=0; 
if v>=0; % upper half space (lower half space zero) 
E=sin(Nz*xsiz/2)/sin(xsiz/2)*sin(Nx*xsix/2)/sin(xsix/2); 
if E==NaN, E=Nx*Nz; end 


end 

if dipdir=='z' 
EFth=2*sqrt(1-w’2)*sin(beta*h*v); % colinear dipoles (z directed) 
EFph=0; 

end 

if dipdir=="x' 


EFth=2*cos(phr).*w*sin(beta*h*v); % parallel case (x directed) 
EFph=-2*sin(phr)*sin(beta*h*v); 

end 

ffth=abs(EFth*E); 

ffph=abs(EFph*E); 

fthit,jt)=ffth; 
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fphdt,jt)=ffph; 
Xthiit, jo=fthatjp)*u; 
Yth(it,jO=fthatjt)*v; 
Zth(it,jO=fthat,jt)*w; 
Xphit,jt)=fphit,jt)*u; 
Yph(it,jt)=fphat,jt)*v; 
Zph(it,jt)=fphdt,jt)*w; 
end 
end 
%plot theta component 
MA X=max([max(max(Xth)),max(max(Yth)),max(max(Zth)),max(max(Xph)),max(max(Yph)),ma 
x(max(Zph))]); 
figure(3),clf 
mesh(Xth/MAX, Yth/MAX,Zth/MA X) 
axis([-1,1,-1,1,-1,1]) 
axis square 
hold on 
% spread N elements along the z axis between [-1,1] 
z0=[-(Nz-1)/2:(Nz-1)/2]*.75/Nz; 
x0=[-(Nx-1)/2:(Nx-1)/2]*.25/Nx; 
[X,Z]=meshgrid(x0,z0); Y=ones(size(X))*.05; 
plot3(X, Y,Z,'ko') 
xlabel('X=IElsin(t)cos(p)’) 
ylabel('Y =IElsin(t)sin(p)') 
zlabel('Z=|Elcos(t)'); 
if dipdir=="x', title(E_\theta for x directed dipoles (parallel case)'), end 
if dipdir=='z', title(E_\theta for z directed dipoles (collinear case)’), end 


%plot phi component 

figure(4),clf 

mesh(Xph/MAX, Yph/MAX,Zph/MAX) 

axis([-1,1,-1,1,-1,1]) 

axis square 

hold on 

% spread N elements along the z axis between [-1,1] 
Z0=[-(Nz-1)/2:(Nz-1)/2]*.75/Nz; 

x0=[-(Nx-1)/2:(Nx-1)/2]*.25/Nx; 

[X,Z]=meshgrid(x0,z0); Y=ones(size(X))*.05; 

plot3(X,Y,Z,'ko') 

xlabel(X=IElsin(t)cos(p)’) 

ylabel('Y =IElsin(t)sin(p)') 

zlabel('Z=|Elcos(t)'); 

if dipdir=="x', title(E_\phi for x directed dipoles (parallel case)'), end 
if dipdir=='z', title(E_\phi for z directed dipoles (collinear case)'), end 
end 
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APPENDIX C. MATLAB CODES FOR BILOTTY?’S METHOD 


This Matlab code was used to compute the voltage ratio versus frequency using the formulas 
derived in Chapter V. 


% FM quadrature demod 
% from note by Coffing and Main 
clear 
rad=pi/180; 
L=.45e-9; Cp=Se-12; Cs=Se-12; Rp=560000; 
cc=Cs/(Cs+Cp); 
wn=1/sqrt(L*(Cs+Cp)); 
Q=Rp/wn/L; 
fn=wn/2/pi; 
disp([{'Cs = 's;num2str(Cs)]) 
disp({'Cp = ',;num2str(Cp)]) 
disp({'L = ',;num2str(L)]) 
disp({'fn = 'num2str(fn/1e6),’ MHz')) 
disp({'Q = ',;num2str(Q)]) 
% frequencies in MHz 
f1=2365; 
f2=2380; 
df=1; 
nf=floor((f2-f1)/df)+1; 
for i=1:nf 
fmhz(i)=f1+(i-1)*df; 
f=fmhz(i)* 1e6; 
w=2*pi*f; % omega 
S=j*w; 
WWw=w/wn; 
Wwp=w+wn; 
wm=w-wn; 
H(i)=j*Q*cc*ww/(1+j*Q/ww*(wp*wm/wn’?2)); 
a=2*Q*wm/wn; 
H1(i)=Cs*L*s%2/((Cp+Cs)*L*s%2+L*s/Rp+1); % equation (2) 
end 
Hm=abs(H); 
Ha=angle(H)/rad; 
H1m=abs(H1); 
H1la=angle(H1)/rad; 


figure(1) 

subplot(121) 
plot(fmhz,Hm) 
xlabel('Frequency, MHz’) 
ylabel(‘lv_2/v_11, dB') 
grid 

subplot(122) 
plot(fmhz,Ha) 
xlabel('Frequency, MHz’) 
ylabel( Angle v_2/v_1, Degrees’) 
grid 


at 


figure(2) 

subplot(121) 
plot(fmhz,H1m) 
xlabel('Frequency, MHz’) 
ylabel(‘lv_2/v_ 11, dB') 
grid 

subplot(122) 

plot(fmhz,H 1a) 
xlabel('Frequency, MHz’) 
ylabel(' Angle v_2/v_1, Degrees’) 
grid 
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